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Motivated by recent Pierre Auger result on the correlation of the highest-energy cosmic 
rays with the nearby active galactic nuclei, we explore possible ultrahigh energy cosmic ray 
^ ■ (UHECR) source distributions and their eflFects on GZK horizons. Effects on GZK horizons by 

^ ' local over-density of UHECR sources are examined carefully with constraints on the degree of 

C/3 local over-density inferred from the measured UHECR spectrum. We include the energy call- 

bration effect on the Pierre Auger data in our studies. We propose possible local over-densities 
of UHECR sources which are testable in the future cosmic ray astronomy. 

^ ■ 1 Introduction 

m 

en ' Recently, Pierre Auger observatory published results on correlation of the highest-energy cosmic 

■ rays with the positions of nearby active galactic nuclei (AGN)^. Such a correlation is confirmed 

. by the data of YakutskEl while it is not found in the analysis by HiResEl. In the Auger result, 

\ the correlation is maximal for the threshold energy of cosmic rays at 5.7 x 10^^ eV, the maximal 

distance of AGN at 71 Mpc and the maximal angular separation of cosmic ray events &tij) = 3.2°. 
Due to increasing efforts on verifying the Auger result, it is worthwhile to examine the above 
correlation from a phenomenological point of view. 
d \ Since the angular scale of the observed correlation is a few degrees, one expects that these 

cosmic ray particle s are predominantly light nuclei. The effect of GZK attenuations on these 
cosmic ray particles'^'^ can be described by a distance scale referred to as "GZK horizon" . By 
definition, the GZK horizon associated with a threshold energy E'th is the radius of a spherical 
region which is centered at the Earth and produce 90% of UHECR events arriving on Earth 
with energies above £'th- 

Assuming a uniform distribution of UHECR sources with identical cosmic ray luminosity 
and spectral index^^, the GZK horizon for protons with E'th = 57 EeV is about 200 Mpc while 
the V-C catalog used by Pierre Auger for the correlation study is complete only up to 100 Mpc. 
Such a deviation may arise from non-uniformities of spatial distribution, intrinsic luminosity and 
spectral index of local AGN as mentioned in^. In addition, the energy calibration also plays 
a crucial role since the GZK horizon is highly sensitive to the threshold energy E^h- Energy 
values corresponding to the dip and the GZK cutof f of UHECR spectrum were used to calibrate 
energy scales of different cosmic ray experiments ^S^. It has been shown that all measured 
UHECR energy spectra can be brought into good agreements by suitably adjusting the energy 
scale of each experiment^. Furthermore, it has been shown that a different shower energy 
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reconstruction method infers a 30% higher UHECR energy than that determined by Auger's 
fluorescence detector-based shower reconstruction^^. 

In this presentation, we report our results^I^ on examining the consistency between Auger's 
UHECR correlation study and its spectrum measurement. The impact by the local over-density 
of UHECR sources is studied. We also study the energy calibration effect on the estimation of 
GZK horizon and the spectrum of UHECR. Certainly a 20% — 30% upward shift on UHECR 
energies reduces the departure of theoretically calculated GZK horizon to the maximum valid 
distance of V-C catalog^. The further implications of this shift will be studied in fittings to the 
shifted Auger spectrum. 



2 GZK horizons and the UHECR spectrum 



Table 1: GZK horizons of UHECR calculated with the local over-density n{l < 30Mpc)/no = 1, 2, 4, and 10, 
and arrival threshold energy Eth = 57 EeV, 70 EeV, 80 EeV and 90 EeV respectively. The listed numbers are in 

units of Mpc. 



n{l < 30Mpc)/no 


Eti, = 57 EeV 


^th = 70 EeV 


Eth = 80 EeV 


Eti, = 90 EeV 


1 


220 


150 


115 


90 


2 


210 


140 


105 


75 


4 


195 


120 


85 


60 


10 


155 


85 


50 


30 



GZK horizons corresponding to different local over-densities and -Eth are summarized in Table 
I. Within the same Eth, local over-densities up to n{l < 30Mpc)/no = 4 do not significantly 
alter GZK horizons. One could consider possibilities for higher local over-densities. However, 
there are no evidences for such over-densities either from astronomical observations^^^l or from 
fittings to the measured UHECR spectrum. We note that GZK horizons are rather sensitive to 
Eth- Table I shows that GZK horizons are ~ 100 Mpc or less for Eth ^ 80 EeV. 

Fittings to the Auger spectrum have been performed in^l^. In our work, we take into account 
the over-density of UHECR sources in the distance scale I < 30 Mpc. The local over-density 
of UHECR sources affects the cosmic-ray spectrum at the highest energy, especially at energies 
higher than 5 • 10^^ eV. Hence the degree of local over-density can be examined through fittings 
to the measured UHECR spectrum. 



Table 2: The values of total from fittings to the Auger measured UHECR spectrum. Numbers in the parenthesis 
are values from fittings to the 8 data points in the energy range 19.05 < logtQ(i?/eV) < 19.75. The last 4 data 

points record events with energy greater than 71 EeV. 



n{l < 30Mpc)/no 


1 


2 


4 


10 


7 = 2.5 


14.12(9.34) 


14.61(9.93) 


17.09(10.50) 


28.09(13.93) 


7 = 2.6 


16.64(12.28) 


15.56(11.90) 


16.01(11.83) 


20.76(11.67) 



Table 3: The total values from fittings to the Auger measured UHECR spectrum with a 30% upward shift on 
UHECR energies. Numbers in the parenthesis are values from fittings to the 8 data points in the energy range 
19.16 < logiQ(i?/eV) < 19.86. The last 4 data points record events with energy greater than 92 EeV. 



n{l < 30Mpc)/no 


1 


2 


4 


10 


7 = 2.4 


8.65(4.30) 


7.39(4.67) 


10.26(6.35) 


27.31(13.34) 


7 = 2.5 


11.82(6.16) 


8.67(5.49) 


7.78(5.23) 


16.18(7.39) 
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Figure 1: Left and right panels depict fittings to tfie original and shifted Auger UHECR spectra respectively with 
the red, green, blue and black curves denote models with local over-densities n{l < 30Mpc)/no = 1, 2, 4, and 10 
respectively. Solid and dash curves correspond to 7 = 2.6 and 7 — 2.5 respectively. We take the source evolution 

parameter m = 3 throughout the calculations. 

The left paenl in Fig. [T] shows our fittings to the Auger measured UHECR spectrum with 
7 = 2.5 and 2.6 respectively. We take the red-shift dependence of the source density as n{z) = 
no(l + z)"^ with m = 3. We have fitted 12 Auger data points beginning at the energy 10^^ eV. 
We make a flux normalization at 10^^ eV while varying the power index 7 and the the degree 
of local over-density, n{l < 30Mpc)/no- Part of values from our fittings are summarized 
in Table II. We found that 7 = 2.5, n{l < 30Mpc)/no = 1 gives the smallest value with 
X^/d.o.f. = 1.57. For the same 7, n{l < 30Mpc)/no = 10 is ruled out at the significance level 
a = 0.001. For 7 = 2.6, n{l < 30Mpc)/no = 10 is ruled out at the significance level a = 0.02. 
We note that, for both 7 = 2.5 and 7 = 2.6, the GZK horizon with n{l < 30Mpc)/no = 10, 
Eti, = 57 EeV, m = 3 and E^ut = 1000 EeV is about 155 Mpc. Since n{l < 30Mpc)/no = 10 is 
clearly disfavored by the spectrum fitting, one expects a GZK horizon significantly larger than 
155 Mpc for Eth = 57 EeV. 

We next perform fittings to the shifted Auger spectrum. The results are shown in the 
right panel in Fig. [1] where the cosmic ray energy is shifted upward by 30%. Part of values 
are summarized in Table III. The smallest x"^ value occurs approximately at 7 = 2.4, n{l < 
30Mpc)/no = 2 with x^/d.o.f = 0.82. One can see that x^ values from current fittings are 
considerably smaller than those from fittings to the unshifted spectrum. Given a significance 
level a = 0.1, it is seen that every local over-density listed in Table III except n{l < 30Mpc)/no = 
10 is consistent with the measured UHECR spectrum. It is intriguing to test such local over- 
densities as will be discussed in the next section. We note that, with a 30% upward shift of 
energies, the cosmic ray events analyzed in Auger's correlation study would have energies higher 
than 74 EeV instead of 57 EeV. The GZK horizon corresponding to £^th = 74 EeV is 120 Mpc 
for n{l < 30Mpc)/no = 2 and 105 Mpc for n{l < 30Mpc)/no = 4. 

We have so far confined our discussions at m = 3. In the literature, m has been taken as 
any number between and 5. It is demonstrated that the effect on UHECR spectrum caused by 
varying m can be compensated by suitably adjusting the power index 7^!^. Since GZK horizons 
are not sensitive to 7 and m, results from the above analysis also hold for other m's. 

3 Discussions and conclusions 

We have discussed the effect of local over-density of UHECR sources on shortening the GZK 
horizon. The result is summarized in Table I. It is seen that such an effect is far from sufficient 
to shorten the GZK horizon at E^^ = 57 EeV to ~ 100 Mpc for a local over-density consistent 
with the measured UHECR spectrum. With a 30% energy shift, each cosmic ray event in 
Auger's correlation study would have an energy above 74 EeV instead of 57 EeV. GZK horizons 



corresponding to Eth = 74 EeV then match well with the maximum valid distance of V-C 
catalog. Fittings to the shifted Auger spectrum indicate a possibility for the local over-density 
of UHECR sources. 

We point out that the local over-density of UHECR sources is testable in the future cosmic 
ray astronomy where directions and distances of UHECR sources can be determined. Table 
IV shows percentages of cosmic ray events that come from sources within 30 Mpc for different 
values of E'th and n{l < 30Mpc)/no. Although these percentages are calculated with 7 = 2.4, 
m = 3 and E'cut = 1000 EeV, they are however not sensitive to these parameters. For E'th = 57 

Table 4: Percentages of cosmic ray events originated from sources within 30 Mpc for different values of Eth and 



n{l < 30Mpc)/no. 



n{l < 30Mpc)/no 


Eti, = 57 EeV 


Eth = 70 EeV 


Et], = 80 EeV 


Eti, = 90 EeV 


1 


0.17 


0.27 


0.36 


0.46 


2 


0.30 


0.43 


0.53 


0.63 


4 


0.46 


0.60 


0.70 


0.77 


10 


0.68 


0.79 


0.85 


0.89 



EeV and n{l < 30Mpc)/no = 1, only 17% of cosmic ray events come from sources within 30 
Mpc. For n(l < 30Mpc)/no = 2 and the same Eth, 30% of cosmic ray events are originated from 
sources in the same region. 

In conclusion, we have shown that the deviation of theoretically calculated GZK horizon 
to the maximum valid distance of V-C catalog can not be resolved by merely introducing the 
local over-density of UHECR sources. On the other hand, if Auger's energy calibration indeed 
underestimates the UHECR energy, such a discrepancy can be reduced. More importantly, 
fittings to the shifted Auger spectrum indicate a possible local over-density of UHECR sources, 
which is testable in the future cosmic ray astronomy. 
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